Topological edge states and quantum Hall effect in the Haldane model 
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We study the topological edge states of the Haldane's graphene model with zigzag/armchair 
lattice edges. The Harper equation for solving the energies of the edge states is derived. The results 
show that there are two edge states in the bulk energy gap, corresponding to the two zero points 
of the Bloch function on the complex-energy Riemann surface. The edge- state energy loops move 
around the hole of the Riemann surface in appropriate system parameter regimes. The quantized 
Hall conductance can be expressed by the winding numbers of the edge states, which reflects the 
topological feature of the Haldane model. 
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The integer quantum Hall effect (IQHE), discovered in 
1980 by Klaus von Klitzing [l], is a striking set of macro- 
scopic quantum phenomena observed in a high mobility 
two-dimensional electron gas (2DEG) in a strong trans- 
verse magnetic field (typically, B ~ 1-30T). Soon after 
the Laughlin's famous gauge invar iance argument and the 
treatment of the IQHE as an adiabatic quantum pump 
Q, it was shortly recognized [3] that the Hall conduc- 
tance (j xy at the plateaus can be understood in terms of 
topological invariants known as Chern numbers [4] , which 
are integrals of the /c-space Berry curvatures of the bulk 
states over the magnetic Brillouin zone. While IQHE 
finds its elegant connection through the adiabatic curva- 
ture with bulk topological invariants, Halperin [5] first 
stressed that the existence of the sample edges, which 
produces the current-carrying localized edge states in the 
Landau energy gap, is essential in the Laughlin's gauge 
invariance argument. Hatsugai further developed a topo- 
logical theory of the edge states [6] , in which topological 
invariants are the winding numbers of the edge states on 
the complex-energy Riemann surface (RS). 

The presence of IQHE fundamentally rely on the 
breaking of the time-reversal symmetry (TRS), which in 
the above mentioned works is brought about by imposing 
an external magnetic field on the electrons. Besides this 
external magnetic field, the TRS also can be broken by a 
variety of the other extrinsic or intrinsic mechanisms. A 
most straightforward way is, like what has been carried 
out in the Aharonov-Bohm (AB) effect, the introduction 
of magnetic flux (instead of magnetic field) to the Bloch 
electrons. In virtue of such a way, by using a graphene 
lattice with the complex hopping matrix elements of the 
next-nearest-neighboring honeycomb sites included, Hal- 
dane first showed that the non-zero Chern numbers and 
thus the IQHE can be realized in zero magnetic field [7[. 
In contrast with the cases with external magnetic field, a 
detailed study of the topological edge states of the Hal- 



dane model is still lacking. This issue is stressed in the 
present paper. Also, our study is motivated by the ob- 
servation that in addition to its importance in charge 
IQHE, the spin-doubled Haldane model, in which TRS 
is recovered, has also played a key role in understanding 
the quantum spin Hall effect (QSHE) [To[ and new phase 
of matter [ill ]. 

Our discussion of the topological edge states for the 
Haldane model begins with deriving the Harper equation 
to describe the wave-function transfer relation between 
two edges in a graphene ribbon. It is found that there 
are two edge states in the bulk energy gap, corresponding 
to two zero points of the Bloch function on the complex- 
energy RS. The edge-state energy loops move around the 
hole of the RS, giving rise to nontrivial winding numbers. 
The quantized Hall conductance can be expressed by the 
winding numbers of the edge states, which reflects the 
topological feature of the Haldane model. 

The graphene lattice is composed of two sublattices 
(denoted by the red and blue dots in Fig. 1). The 
graphene ribbons with zigzag edges and armchair edges 
are plotted in Fig. 1(a) and 1(b), respectively. The lat- 
tice tight-binding Hamiltonian [7] is given by 



H 



E** 



(i) 



In the above Hamiltonian, the on-site energy to=+M on 
A site and — M on B site. t\ and ti are real hopping 
matrix elements between nearest neighbors on the dif- 
ferent and the same sublattices, respectively. To break 
TRS, a complex phase <pij is introduced to the next near- 
est neighbor hopping ^. Following Haldane 0], we set 



the magnitude of this complex phase as \(j>i 



and 
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the direction of the positive phase is anticlockwise. Note 
that the net flux is zero in one unit cell. Since the spin- 
orbital effect is not included, we neglect the spin indices 
for simplicity. 

Now, let us derive the Harper equation [H, Q of the 
graphene ribbons with zigzag edges. We suppose that 
the system is periodic in the x direction while it has two 



2 




and ^nA was replaced by ip n . Eq. (|4|) is the so-called 
Harper equation [8|, l2[ . The next key step is to represent 
Eq. (|4]) in the transfer matrix form. After a tedious but 
straightforward derivation, we find that by introduction 
of a new wave function y? n , which is a linear transforma- 
tion of the original wave function 
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<£n = VVi H Wn-l + 1pn-2, (6) 



FIG. 1: (Color online) The structure of graphene ribbons with 
(a) zigzag edges and (b) armchair edges. The rectangle with 
the dashed line is the unit cell. 



edges in the y direction [see Fig. 1(a)]. In the following 
we replace index i with (mns) to denote the lattice site, 
where (mn) label the unit cells and s label the sites A 
and B in this cell. The distance of the nearest neighbor- 
ing lattice sites is set to be unity throughout this paper. 
Since the lattice is periodic in the x direction, we can use 
a momentum representation of the electron operator 
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where (X mns , Y mns ) represents the coordinate of the 
site s in the unit cell (mn), and k x is the momentum 
along the x direction. Let us consider the one-particle 
state \ip(k x ))=^2^ ns (k x )cl ls (k x ) |0). Inserting it into the 
Schrodinger equation H \ip)=e one can easily get the 
following two eigenvalue equations for sites A and B: 

€t/j nA = (M + /+) tpriA + g- W>(n+1)A + ^(u-I)a) 

+ *l^(n-l)£ + 90^nB: 
t^nB = (~M + f-)^nB + 0+ W>(n+1)£ + ^(u-I)b) 

+ tl1p(n+l)A + gO^nA, (3) 

where f±=2t2 cos (V3k x ± (f) , g±=2t2 cos (^^k x ± (j^J , 

and go=2ti cos (^^-k^j. Eliminating the B sites, we ob- 
tain a difference equation for A sites 

VVi = &W>(n-l) + ^(n-3)) + #(n-2) ~ ^(n-4), (4) 

where 
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then the new wave function cp n can be written in the 
following transfer matrix form 
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t -1 
1 



<Pn-l 



M(e) 
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^ (7) 
where ^— b+ V b + 4 ( 2 +0 ^ ]y[ ore generally, we take y?o an d 
(^l as the wave functions at two open edges. Then we 
get a reduced transfer matrix linking the two edges as 
follows: 



"PLy 



M(e) 



(8) 



where 



M (e) = M(e) 



Mn(e) M 12 (e) 
M 2 i(e) M 22 (e) 



(9) 



All kind of solutions of Eq. (|SD are obtained by different 
choices of and 

Similarly, one can obtain the eigenvalue equations of a 
graphene ribbon with armchair edges, 
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+ fg [e-^V(n+i)B + e'*V(n-l)B] • (10) 



However, because the derivation of the Harper equation 
for a graphene ribbon with armchair edges is too sophis- 
ticated, here we do not write out the transfer-matrix ex- 
pression of this Harper equation. Moreover, because the 
main results and the discussions on the graphene ribbons 
with zigzag and armchair edges are similar, in the follow- 
ing we focus our attention to the graphene with zigzag 
edges. The general open boundary condition is 



(ll) 
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With Eqs. ([8]) and ([9]), one can easily get that the solu- 
tions satisfy 



and 



M 21 (e) =0 



<PL v -i = ~Mu (e)cpi. 



(12) 



(13) 



If we use a usual normalized wave function, the state is 
localized at the edges as 

\Mn (e) | <C 1 localized at y « 1 (down edge), 
I Mil ( e ) | ^> 1 localized at ~ L y — 1 (up edge). 

(14) 

Because the analytical derivation of the energy spec- 
trum in the presence of edges is very difficult, we now 
start a numerical calculation from Eqs. (J3]) and ([T0|) . 
Varying all the controllable parameters, which are the 
relative site energy M/ti, the next nearest neighbor hop- 
ping ^Ai, and the complex phase 0, we can get three 
different cases happening in the energy spectrum of the 
graphene ribbons. We draw in Figs. 2(a)-(c) [Figs. 3(a)- 
(c)] the energy spectrum of graphene ribbons with zigzag 
(armchair) edges as a function of k x for these three dif- 
ferent cases, i.e., the case Mjti < 3 a/3 sin (case I), 
the case Mjti < — S^sm (j) (case II), and the case 
Mjti > 3a/3| sin 0| (case III), respectively. The num- 
ber of sites A (B) in y direction is chosen to be L y =40. 
Clearly, from Figs. 2 and 3 one can see that there are 
two dispersed energy bands (the shaded areas) with two 
edge states (the colored lines) lying in the energy gap. 
It is our task to show that the geometric nature of the 
edge states in these three kinds of parameter regions are 
totally different, which can be described by the wind- 
ing numbers of the edge states on a complex energy RS 
within the topological edge theory [6[ . 

To show this, first, we ignore the open boundary con- 
dition and consider the bulk Bloch function at sites with 
^/-coordinate of (L y — 1). For Bloch function, cp^ and 

(p^ compose an eigenvector of M with the eigenvalue p, 



M(e) 




(15) 



In order to discuss the wave function of the edge state, we 
extend the energy to a complex energy. In the following, 
we use a complex variable z instead of real energy e. From 
Eq. $15D we get 



p(z) = \ A(z) 



(16) 



and 
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-M 21 (z) Vu (17) 




FIG. 2: (Color online) Energy spectrum of the graphene rib- 
bon with zigzag edges under different complex phases parame- 
ters: (a) 0=7r/3, (b) (j)—— 7r/3, (c) 0=7t/6. The other param- 
eters are set as M/£i=l and £2/£i=l/3 in all three figures. 
The shaded areas are the energy bands and the colored lines 
are the spectrum of the edge states. The red (solid) and blue 
(dashed) lines mean that the edge states are localized near 
the down and up edges, respectively. 




FIG. 3: (Color online) Energy spectrum of the graphene rib- 
bon with armchair edges under different complex phases pa- 
rameters: (a) 0=7r/3, (b) 0=— 7r/3, (c) </>=7r/6. The other 
parameters are set as M/t\=l and £2/£i=l/3 in all three fig- 
ures. The shaded areas are the energy bands and the colored 
lines are the spectrum of the edge states. The red (solid) and 
blue (dashed) lines mean that the edge states are localized 
near the down and up edges, respectively. 



where A(z)=Ti[M(z)] and cj= v / A 2 (z) - 4. Clearly, 



detM(e) 



(18) 



since detM(e)=l. From Eq. (|TT|) one can find that 
the analytic structure of the wave f unction is deter- 
mined by the algebraic function uu=^/ A 2 (z) — 4. The 
RS of uj—\J 'A 2 (z) — 4 on the complex energy plane can 
be built by the conglutination between different analytic 
brunches. Here, the close complex energy plane can be 
obtained from the open complex energy plane through 
spheral pole mapping [see Fig. 4(a)]. Now let us discuss 
the analytic structure of uj = y A 2 (z) — 4 on the open 
complex energy plane. If the system has q energy bands, 
i.e., 
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with q secants. Whereas, if we set arg(z — A2j-i)=27r 
and arg(z — A2j)=7r, which corresponds to co(jiij-i) < 
(oo(/^j) > 0) when z moves along c\ (02), then the branch 
R~ is defined as a complex plane with q secants. The 
definitions of R + and R~ for odd values of j are reverse 
to those for even values of j. So, if z lies in the jth gap 
from below on the real axis, 

a uo > 0, z: real on R a (a = +, -), 

and at energies fij (eR a , a=+, — ) of the edge states we 
have 

w(jij) = a |M„ - M 22 ( M ,)| , (21) 

In addition, one can easily obtain 



FIG. 4: (Color online) (a) The open complex energy plane are 
mapped to the close complex energy plane through spheral 
pole projection, (b) Two sheets with two cuts which corre- 
spond to the energy bands of the graphene nanoribbons. The 
RS of the Bloch function is obtained by gluing the two spheres 
along the arrows near the cuts. 



e e [Ai, A 2 ] , [A 2 j-i, A 2 j] , [A 2(? -i, X 2q ] , (19) 

where Xj denote energies of the band edges and A^ < Xj , 
i < j, then uo can be factorized as 



2q 



Ai)- 



(20) 



In the present Haldane model, there are two energy 
bands, so q=2. The two single- value analytic branches 
are defined on the same complex energy plane with q 
secants. The difference between the two branches are 
specified in the following paragraph. 

For an up- or down-edge-state energy fij in the gap 
[A2j, A2J+1], In order to ensure uo(jij) > 0, we divide the 
two single-value analytic branches in terms of the par- 
ity (evenness or oddness) of j. Let us consider the case 
that the energy z lies in the band [\2j-i,^2j] [ see Fig. 
4(a)]. The left side of this energy band is the (j — l)th 
gap, while the right side is the jth gap. When the en- 
ergy z moves from the jth band to the (j — l)th gap (the 
jth gap) along an arbitrary path c\ (02), only the sin- 
gularities A2J-1 and A2j have contributions to the vari- 
ance of the principal value of the argument of uo. On 
the up bank of the secant, we distinguish two branches 
R + and R~ as the following: For even values of j, if we 
set arg(z — A2j-i)=0 and arg(z — A2j)=7r, which corre- 
sponds to uo(jiij-i) > {uj(iJLj) < 0) when z moves along 
ci (C2), then the branch i? + is defined as a complex plane 



x f < — 2 for j odd 
1 > 2 for j even ' 



(22) 



where the energy e (on R a ) is in the jth gap. 

When the branches i? + and R~ on the open com- 
plex energy plane are mapped to the close complex en- 
ergy plane through spheral-pole-projection, one can get 
two single- value analytic spherical surfaces. The RS is 
obtained by gluing the two spherical surfaces at these 
branch cuts with making sure that the ± banks face 
the =p banks of other sphere [see Figs. 4(b)]. Note that 
there are two real axes after gluing. In the present model 
the genus of the RS is g=l, which is the number of en- 
ergy gaps. The wave function is defined on the g=l RS 
Y, g= i(k x ). The branch of the Bloch function is specified 
as cj>0, which we have discussed above. With Eqs. 



(|2T|) and (QZj), and using the fact that ^L y -i(/ij)=0 for 
Hj^R a and ipL y -i(lJ>j)¥ z Q f° r l^j^R~ a ^ one can obtain 
that when the zero point is on the upper sheet of RS 
(i? + ), the edge state is localized at the down edge; when 
the zero point is on the lower sheet of RS (R~), the edge 
state is localized at the up edge. 

Figure 5(a)-(c) schematically show the RS's for the 
present Haldane model with the system parameters be- 
longing to case-I (M/t2 < SV^sine^), case-II (M/t2 < 
— 3 a/3 sin cV), and case-Ill (M/t2 > 3\/3| sin0|), respec- 
tively. On each RS T lg= i{k x ) the energy gap corresponds 
to the loop around the hole of the T lg= i(k x ) and the en- 
ergy bands correspond to the closed paths vertical to 
the energy gap loop on T, g= i(k x ). The Bloch function 
is defined on this surface. For the fixed k x and 0, there 
is always g=l zero point at the down-edge-state energy 



(down) 



Since there are two real axes on the T, g= i(k x ), 
correspondingly, there is g=l zero point at the up-edge- 
state energy fi^^ . 

The above analysis is for the fixed k x . Changing k x in 
one period, we can consider a family of RS's T, g= i(k x ). 
T ig= i(k x ) can be modified by this change. However, all 
the RS's T, g= i(k x ) with different k x are topologically 
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FIG. 5: Riemann surfaces of the Bloch functions for different 
winding numbers: (a) 1=1, (b) I=— 1, and (c) 1=0. 



equivalent if, as what happens in the present model, 
there are stable energy gaps in the 2D energy spec- 
trum. By identifying the topologically equivalent RS's 
^g=i{k x ), the behavior of the track of Hj(k x ) (including 

the up-edge-state energy ji^^ and down-edge-state en- 
ergy /i^ down ^) depends on system parameters, as shown 
in Fig. 5. In Fig. 5(a), which corresponds to case-I 
of M/t 2 < 3\/3sin0, one can observe that by varying 
k X: the down-edge-state energy /i^ down ^ (k x ) moves from 



the lower band (band 1 in Fig. 5) edge to the upper 
band (band 2 in Fig. 5) edge, while the down-edge- 
state energy /i^ down ^ (k x ) moves from the upper band edge 
to the lower band edge. That is to say, the two edge- 
state energy tracks in the same energy gap move around 
the hole and form an oriented loop C(jHj). In case-II of 
M/t2< — 3-\/3sin(/>, as shown in Fig. 5(b), the two edge- 
state energy tracks moving around the hole also form 
an oriented loop. However, the orientation of the loop 
is right about with respect to that in case-I. Finally, in 
case-Ill of M/t 2 >3\/3\ sin</>|, as shown in Fig. 5(c), one 
can observe that the fij(k x ) moves along the hole and 
turns back before arriving at the second energy band. In 
this case, the two edge-state energy tracks in the same 
energy gap approximately form two circularities. 
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I(A,B)=+1 I(A,B)=-1 
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FIG. 6: Intersection number I(A,B) of two curves A and B. 
Each intersection point contributes by +1 or —1 according to 
the direction. 



It is known that on a general genus-g RS, all kinds 
of loops (the first homotopy group) are generated by 2g 
canonical loops (generators), oti and ft, z=l, . . . ,g. See 
Fig. 5 for g=l. The intersection number of these curves 
(including directions) [6[ is given by (see Fig. 6) 



Sij. 



(23) 



Any curves on the RS are spanned homotopically by oti 
and Pi. When the edge-state energy loop fij(k x ) moves 
p times around the jth hole with some integer p, one has 



which means 



I(a.i,C(jij)) =p5 k 



(24) 



(25) 



When the Fermi energy ep of the 2D system lies in the 
jth bulk energy gap, the Hall conductance is given by 
the winding number of the edge state [6], which is given 
by the number of intersections I(aj, C(fij)) (= I(C(fij))) 
between the canonical loop ctj on the RS and the trace of 
H j. In the present single-gap model, we obtain the Hall 
conductance provided by the edge states as follows 



edge 



h 



■J(CM). 



(26) 



From Figs. 5(a) and (b) it can be observed that /i moves 
one time across the hole (p=l), which in terms of Eq. 
d23j) means C(fi) « (3. Considering simultaneously the 
winding direction (see Fig. 6), one can obtain that in 
Fig. 5(a) 7(C(/i))=-l, while in Fig. 5(b) J(C(/i))=l. 
In Fig. 5(c), because p=0, / (c^, C(/ii))=0 and the Hall 
conductivity is zero. Therefore, we get 



edge 



h ' 
h ' 



M/t 2 < 3V3 sin 6 
M/t 2 > 3V3\ sin 0| 
M/t 2 < -3V3 sin 6 



(28) 



Finally, let us compare this result for the graphene rib- 
bons with zigzag/armchair edges with the bulk graphene, 
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in which the topological invariant is the Chern number. 
In the bulk Haldane model, when the Fermi energy ep 
lies in the energy gap, the Hall conductance is quantized 
as cr^y lk = \C\, where C\ is the Chern number of the 
lower energy band. It turns out that 

{1, M/t 2 < 3V3 sin 6 
0, M/t 2 >3V / 3|sin0| . (29) 
-1, M/t 2 < -3v/3sin^ 

From Eqs. ([28]) and ((29]) one can obtain a^ ge = o^ lk , 
which is in accord with the established recognition [H[ Q 
on the Hall conductance in the systems with and without 



edges. 

In summary, we have investigated the topological prop- 
erty of the edge states in the Haldane model. The Harper 
equations for solving and analyzing the edge states have 
been derived. It has been found that there are two edge 
states lying in the bulk energy gap. These two edge states 
move with varying k x around the hole in the RS and form 
an orientated energy loop. With the winding number 
of the edge states, we have obtained that the edge- state 
Hall conductance is cr|^ ge =±-^- or under different cases, 
which agrees with that based on the topological bulk the- 
ory. 



[1] K. von Klitzing, G. Dorda, and M. Pepper, Phys. Rev. 

Lett. 45, 494 (1980). 
[2] R. B. Laughlin, Phys. Rev. B 23, 5632 (1981). 
[3] D. Thouless, M. Kohmoto, M. Nightingale, and M. den 

Nijs, Phys. Rev. Lett. 49, 405 (1982). 
[4] D.J. Thouless, Topological Quantum Numbers in Nonre- 

latisvistic Physics (World Scientific, Singapore, 1998). 
[5] B. I. Halperin, Phys. Rev. B 25, 2185 (1982). 
[6] Y. Hatsugai, Phys. Rev. Lett. 71, 3697 (1993); Phys. 

Rev. B 48, 11851 (1993). 



[7] F. D. M. Haldane, Phys. Rev. Lett. 61, 2015 (1988). 
[8] P. G. Harper, Proc. Phys. Soc. London Sect. A 68, 874 
(1955). 

[9] D. R. Hosfstadter, Phys. Rev. B 14, 2239 (1976). 
[10] C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 226801 
(2005). 

[11] C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 146802 
(2005). 



